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Abstract 

Oscillatory tunneling magnetoresistance (TMR) as a function of spacer thickness is investigated 
theoretically for a magnetic tunnel junction with a nonmagnetic layer inserted between the tunnel 
barrier and the ferromagnetic layer. TMR is characterized in an analytical form, that is expressed 
with the transmission and reflection amplitudes of single interfaces at the Fermi level, and by the 
extremal wave vectors. Electronic structures with multiple bands are taken into account in the 
derivation characterizing the TMR, and the proposed analytical expression can be directly ap- 
plied to real junctions. Based on our model, the features of TMR dependence on spacer thickness 
are discussed, including selection rules for the oscillation period. Numerical calculations are per- 
formed using an envelope-function theory for several cases, and we show that our model is in good 
agreement with the exact result. 

PACS numbers: 72.25.Hg, 72.25.Mk 
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I. INTRODUCTION 



Since high tunneling magnetoresistance (TMR) was first observed at room temperature,- 
magnetic tunnel junctions (MTJs) have been a focus of interest. Extensive research has 
been carried out to understand and improve the properties of MTJs. A huge increase in 
TMR with lower junction resistance was achieved when A10 x . tunnel barriers were replaced 
by MgO, and this was followed by realizations of memory devices based on MTJs with MgO 
barriers.- Tunneling current in the MTJ is spin polarized, which adds another dimension 
to the tunneling effect, and scientific attention has thus been drawn to the spin-dependent 
tunneling phenomenon. When a nonmagnetic (NM) layer is inserted between a ferromagnetic 
(FM) layer and the insulating (I) tunnel barrier of the MTJ, the spin polarization of the 
tunneling current changes and the TMR is directly affected. An early theoretical work 
predicted the oscillatory TMR as a function of the NM thickness due to the quantum well 
states inside the NM layer. 3 In sputtered samples, it has been shown that an NM layer 
between the tunnel barrier and FM layer could be detrimental to TMR, and the TMR 
decreases as a function of NM thickness.- 1 ^ These experimental results have been explained 
theoretically with a free electron model, and the decay of TMR was attributed to a loss of 
coherence in the electron propagation. 6 Different experimental results have been obtained for 
a crystalline NM layer inserted between the tunnel barrier and the FM layer. Yuasa et ai— 
experimentally investigated the dependence of TMR on Cu thickness in NiFe/A10 x /Cu/Co 
junctions with samples grown by molecular beam epitaxy. They found that the TMR decayed 
but oscillated as a function of Cu thickness. The oscillation period was determined by 
the nesting feature of the Cu Fermi surface. This oscillatory TMR has been investigated 
theoretically based on a single-band tight-binding model, a free-electron model, and full- 
band calculations.- - — Many features have been explained with calculations using simple 
models, but direct comparisons with the experimental data are difficult because realistic 
electronic structures were not considered. The full-band calculations are very useful for the 
description of real systems. However they are time consuming, and sometimes it is not 
easy to understand the underlying physics. Furthermore, full-band calculations are usually 
carried out for an ideal situation, and significant discrepancies often occur between theory 
and experiments. 

We introduce an analytical expression that describes the dependence of TMR on NM 
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thickness for FM/I/NM/FM junctions based on full-band structures. Our approach uses 
the generalization of a previously described single-band case^ 9 to a multiple-band case that 
considers the real materials. The TMR is expressed with transmission and reflection am- 
plitudes of single interfaces at the Fermi level, and extremal wave vectors. The full-band 
structures of the materials are taken into account in our proposed model, and the calcula- 
tion of several transmission and reflection amplitudes with real band structures can make a 
direct comparison with experimental results possible. Based on our model, selection rules 
for the oscillation period are discussed, and we suggest that very few oscillation periods 
will be observed in experiments even when there are many extremal spanning vectors of the 
NM Fermi surface. This situation is very different than the interlayer exchange coupling in 
magnetic multilayers. Our model explicitly shows that TMR dependence on NM thickness 
is affected by the thickness of the tunnel barrier, and predicts that the TMR will go to zero 
as the NM thickness increases. To check the validity of our model, we carried out numerical 
calculations using an envelope-function theory for several cases including NM material with 
a Fermi surface similar to the Cu(OOl) case. Although our model is calculated using the 
parameters at the Fermi level and the extremal wave vector, it is in good agreement with 
the exact numerical results. 



II. THEORETICAL MODEL 



Figure 1(a) shows a schematic diagram of an MTJ with an NM layer inserted between the 
tunnel barrier and the right magnetic layer [FM(R)]. The growth direction is taken as the 
z axis, and d is the thickness of the NM layer. We used a frozen potential approximation, 
and the eigenstate of the MTJ is expressed with linear combinations of the bulk states for 
each layer. We assumed that the wave vector component parallel to the interface (ky) is 
conserved throughout the MTJ. A two-channel model was adopted and spin-flip scattering 
was ignored. The normalized bulk solution of the material for the left magnetic layer [FM(L)] 
is denoted as e, ky, kz%a for a given energy e and ky, where k z is the z component of the 
wave vector, n is the band index, a is the spin index, and the + (— ) sign is for the state 
traveling to the right (left). Similarly, the bulk solutions of the FM(R) and NM materials 



are expressed as 



e, ky, kf^ ) ) and 



e, l&H,k^n ^ ), respectively. Multiple bands are taken 



into account, and 2Nl(r)o- is the number of bulk states in the left (right) FM layer for a 
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given e, kn, and spin a. The number of bulk states in the NM layer for a given e and kii is 
denoted as 2N^- The eigenstate of the MTJ is written as 
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where v z is the z-component of the group velocity [v z = (l/h)(de/dk z )] for the corresponding 
bulk eigenstate, and At, Ct, and B^ a are coefficients to be determined from the boundary 
conditions. Note that the bases in Eq. ([I]) are adjusted so that the current is normalized. 
The eigenstate inside the tunnel barrier is not shown here. We define vectors and 
as A± = (At, A± L J T and B± ee (£?± , £f± , ■ • ■ , B± R J T . Then, A; and B+ are 

related to A+ and B^ by the S'-matrix,— 



(2) 



Matrix element t a ^ nn i (r a ^ nn i) is a kind of transmission (reflection) amplitude for an incoming 
wave from the left \e, ky, k^, a ) to be transmitted (reflected) to \e, ky, kft) (\ £ , ky, k^~ a ~)). 
The transmission (reflection) amplitude for the opposite direction is given by t' ann , (r' ann ,). 
We calculated the conductance G for low bias and zero temperature from the Landauer- 
Biittiker formalism as follows 

G=^^Tr[tt(£ F ,k||)t (T (e F ,k||)], (3) 





where Bp is the Fermi energy. As shown in Figs. 1(b) and (c), we considered the FM(L) /I/NM 
and NM/FM(R) interfaces separately, and expressed t a of the MTJ with the transmission 
and reflection amplitudes of each separated interface. The eigenstate of the FM(L)/I/NM 
system shown in Fig. 1(b) is expressed as 
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where and C^f are coefficients. We define = (A£, Ag 1 ,---, J t and C£± 



(Cf^, C^, ■ • • , Cjv^cr) T ' an d the relation between the coefficients is expressed as 
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Similarly, for the NM/FM interface shown in Fig. 1(c), the eigenstate is given by 
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The vectors C R± = (C*±, C*±, • • • , C*$ a ) T and = (Sf±, B*± • • • , fljg J T are related 



as follows 
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r and t of the MTJ in Eq. (El) can be expressed with r L , t L , r' L , and t' L in Eq. (EJ and r^, 
t H , r ,i? , and t' R in Eq. (J7J) by considering the multiple reflection inside the NM. We introduce 
the mean free path A due to scattering inside the NM layer. For simplicity, we assumed that 
A is constant, although the dependence of A on other parameters can be included in our 
calculation. Then, the phase- coherent part of the reflection amplitude r c is given by 

oo 

r c = r L a + t'^e-x + £ ^{p R J L a Y p R r^^\ (8) 



n=l 



where the matrix elements p R nn „ r^ nn ,, and p' a ^ nn , are p R nn , = e lk ^ d r R nn ,, r^ nn , = 
e ik z +d t^ nn , } and p'^nn' = e lhz '" d r'^ nn ,. The phase-coherent part of the transmission ampli- 
tude t c can be obtained in a similar way. Because of the scattering inside the NM layer, we 
have Tr[t c H c ] + Tr[r c ^r c ] < 1 and we need to include the diffusive part of the transport. We 
assume that the transmission back to the FM(L) layer through the tunnel barrier is much 
smaller than that through the NM/FM(R) interface. Then, the 1 - Tr[t cf t c ] - Tr[r ct r c ] 
portion contributes to the sequential transmission.- 2 Adding the coherent and sequential 
transmissions, we have Tr[tH] = 1 — Tr[r c ^r c ]. Finally, using the properties of the S matrix 
and taking the first-order term in e~ 2d / A , we obtain 

Tr[t f t] = Tr[t Lt t L ] + 2Re Tr^V^T^e-x. (9) 



Even when A is very large, this is a reasonable approximation because the magnitude of the 
matrix element r^ nn , is less than 1, and the higher-order terms are more rapidly oscillating 
as functions of d and consequently contribute less to the conduction. The conductance is 
given by 

G = G + ^e-xRe£Tr[t^V^], (10) 



h 
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where Go is the conductance of the FM(L)/I/NM junctions and the energy is set to the 
Fermi level (e = ep). The conductance depends on the magnetic configurations, and we 
denote the conductance for parallel (anti-parallel) magnetization of two magnetic layers 
as (jp(AP)- The TMR is given by AG/Gap, where AG is AG = Gp — Gap- Here, we 
will show the calculation of AG, and Gp(ap) — Go can be obtained in the same way. We 
define AT£, and Ar* , as AT*, = (tftfVf - tftfVf) = |AT£,| e^«' an d Ar* , = 

\ / ran' 

{ r fnri ~ r inn') /2 = \ Ar^ n , \ ^ nn ' , where | (|) is the majority (minority) spin. AG is 
expressed as 

AG = ye-TRe^^ \AT^\ |Ar* ,| e^'^w), (n) 

ra, n' k||,(T 

where g nn / and <j) nn , are g nn / = k% + - k%~ and </w = <p^ n , + 0^ n ,. 

The summation over kn can be performed in a manner similar to the calculation of the 
interlayer exchange coupling in magnetic multilayers.— 1 ^ |AT^ n , | and e tq ™' d are rapidly 
changing as functions of ky. We assume the exponential dependence of lAT^, such 
that lAT^J cx e~ bXnn '^\ Suppose that (k ax ,k ay ) is an extremal point, which means 
Vit, [—bxnri + i{<inn'd + 4>nn')] = at k|| = (k ax ,k ay ). Since the main contribution to the 
integral comes from the vicinity of the extremal point, we expand —bxnri + i(<lnn'd + </w) 
around the extremal point as follows 

- bXnn' + i(qnn'd + 4> nn >) ~ ~b Xa + * + 4>a) 

Zb % ~ Zd )( k x-k ax ) 2 (12) 
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where new parameters -jj- = | , -5- = I ^ila"' , -3- = I ^It' , 4- = I ^It' , d QX = 
\ K ax d dk 2 "' anc ^ day = \ K a y & tk 2 are evaluated at e = Bp and k|| = (k ax ,k ay ). Then, the 
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summation over ky is carried out analytically and AG is given by 

p 2 n I AT L \ \ Ar R \ pKiad+^a) 
till ^ b • d+d ax b • d+d a y 

V ^ax ^ax y ^ay ^ay 

where n a is the number of the extremal points of the same kind. The phase of the square root 
is taken from —tt/2 to n/2. The parameters in Eq. (Tl3l) are evaluated at the Fermi level and 
the extremal point. Suppose we have Vky X — at kjj = (k x , k b y ) and Vky (qd+4>) = at ky = 
(k x , ky). In general, kjj is different from kjj and the corresponding extremal point (k ax , k ay ) 
is a complex number. This makes other parameters such as q a complex numbers, and the 
situation is rather complicated. However, when kjj and kjj are far apart, the contribution 
is negligible because of small |AT^|. The most important case is when kjj and kjj coincide. 
This is expected to happen often at ky = due to symmetry. In this case, k ax and k ay are 
real and AG becomes 
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where we used rj L = r| L with |r| L | = 1, which is expected for a typical thickness of the 
tunnel barrier. In this case, q a is exactly the extremal spanning vector of the NM Fermi 
surface. 

The extremal spanning vector q a of the NM Fermi surface gives rise to the period of the 
TMR oscillation. In principle, multiple periods are possible depending on the shape of the 
NM Fermi surface. However, the extremal spanning vectors of the NM Fermi surface are 
the necessary condition for TMR oscillation, and the period would be observed only when 
the corresponding |AT^| is large enough. Except for the symmetry point in ky space, the 
extremal points of the NM Fermi surface hardly coincide with the maxima point of |AT^| 
of the tunnel barrier. Thus, only few periods would be observable among the possible 
periods from the extremal spanning vectors. Compared with the oscillation of the interlayer 
exchange coupling for a corresponding spacer, fewer periods would be observed in TMR 
oscillation. Moreover, as b increases, |AT^| decreases rapidly as a function of ky away from 
the maxima. Thus, the multiple periods would be more difficult to observe for a thicker 
tunnel barrier. Since AT^I is material dependent, the period of TMR oscillation could 
change for a different tunneling barrier. In Eq. (Q3J), the extra phase factor arises from the 
square root in the denominator. This extra phase factor is d dependent, and the measured 



oscillation period will be slightly different from 2ir/q a . For instance, the period of TMR 
oscillation for NiFe/AlO x /Cu(001)/Co is expected to be slightly smaller than 27r/q a , which 
is the period of the interlayer exchange coupling for fee Co/Cu(001)/Co multilayers. 

We considered the observable periods with the assumption that the tunneling probability 
is maximum at kn = 0. Then, only the oscillation for the extremal spanning vector at 
k|| = would be observed. In NiFe/AlO :r /Cu(001)/Co junctions, there are two possible 
periods from the Fermi surface of Cu along the (001) direction. The long period is from 
the extremal spanning vector at kn = , and the extremal vector for the short period is far 
away from kn = 0. Thus, the short period is intrinsically invisible, and only the long period 
could be observed. There are four possible oscillation periods in NiFe/AlO :r /Cu(110)/Co 
junctions from the Cu Fermi surface along the (110) direction.— Among them, the only 
observable period is the extremal spanning vector corresponding to ky = 0. This period is 
very short and would be easily wiped out by the interface roughness. Thus, it would be 
difficult to observe a TMR oscillation in this system. There is only one possible period for 
Cu(lll)^ i 2^ and this period has been observed experimentally in the interlayer exchange 
coupling of Co/Cu(lll)/Co multilayers. 2 - However, this period would not be observed in 
NiFe/AlO-r/ Cu(lll) / Co junctions because the corresponding extremal spanning vector is far 
from ky = 0. Therefore, except for the long period of NiFe/AlO^/Cu(001)/Co junctions, it 
would be difficult to observe a TMR oscillation as a function of Cu thickness. The situation 
for the Au and Ag spacers will be same as in the case of Cu because the shapes of the Fermi 
surfaces are similar. In the Fe/Cr/Fe multilayers, the long periods of interlayer coupling 
oscillation have been observed as a function of the Cr thickness for the (100), (110), and 
(211) orientations.— However, these long periods would not be observed experimentally in 
Fe/A10 x /Cr/Fe or Fe/MgO/Cr/Fe tunnel junctions because the long periods are from the 
iV point of the Cr Fermi surface,— 1 ^ which is away from kn = for any orientation. In 
experiments, this long period has not been observed for Fe/AlO x /Cr(001)/Fe (Ref. l26l ) and 
Fe/MgO/Cr(001)/Fe (Ref. |27|) tunnel junctions. Usually, the long period is clearly observed 
because it is not eliminated by the interface roughness. Although we discussed the case that 
the tunneling is dominated by perpendicularly incident electrons, the analysis is similar when 
the tunneling probability is maximum or high away from kn =0. If the dependence of the 
tunneling probability on kn changes with a different tunnel barrier, the TMR dependence 
on the NM spacer will be altered accordingly. Still, the crucial criterion is whether the 
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tunneling probability is significant or not at the point of the extremal spanning vectors of 
the NM Fermi surface in the ky space. Except special points, the chances are that the point 
of the maximum tunneling probability does not coincide with the position of the extremal 
spanning vectors of the NM Fermi surface in the ky space. Thus many oscillations inferred 
from the NM Fermi surface would not be observed. When the tunnel barrier is extremely 
thin, the situation can be much different because the tunneling probability dependence on 
ky may change significantly. More oscillation periods can be observed with thinner tunnel 
barriers. Even in this case, the oscillation periods associated with relatively high tunneling 
probability will be observed. 

Without the scattering effect, AG and TMR decays as 1/d for a thick NM layer. How- 
ever, for thin NM layers, the decay rate is slower than 1/d and is affected by the tunnel 
barrier thickness b, and also by d ax and d ay (ky dependence of the reflection-amplitude phase 
factors). The amplitude decays much more slowly than 1/d for coherent transport when d 
is of the same order of magnitude as b. In experiments, the TMR oscillation decays much 
faster than 1/d, which seems to be due to scattering. As the NM thickness increases, our 
model predicts that AG and TMR go to zero even when the mean free path (A) is very long. 
We will address this point in Sec. IIIC. 

III. NUMERICAL CALCULATION WITH AN ENVELOPE-FUNCTION THE- 
ORY 

To test the validity of our model, we carried out numerical calculations based on an 
envelope-function theory for several cases. We used the same material for FM(L) and FM(R), 
and ignored scattering. The continuity of the wave function and the conservation of current 
at the interface were taken as the boundary conditions. 

A. Effective- mass approximation 

First, we considered the case that the dispersion relation of the NM material is the same 
as that of the FM material for majority spin. The dispersion relations of the FM material 
and the insulator are given by e(k) = h 2 k 2 /2m* a + V a and e(k) = h 2 k 2 /2m Q + Vi, respectively, 
where m* is the spin-dependent effective mass of the FM material, m is the bare electron 
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mass, and Vi is the height of the tunnel barrier. We set = [Vj. = A] for the majority 
(minority) spin in the FM layer using the spin-splitting energy A. The effective mass of the 
FM material is m-j- = m for the majority spin and ml = ttlqE f / (sp — A) for the minority 
spin. Schematics for the dispersion relations of the FM, I, and NM materials are shown in 
Fig. 2. Under these conditions, all the traveling states in the NM layer have corresponding 
traveling states in the FM layer, and total reflection does not occur at the NM/FM interface 
for any kn. The parameters used in the calculation are Ef = 4 eV, Vj = 6 eV, and A = 2.5 
eV. The TMR is plotted as a function of NM thickness d for the tunnel-barrier thicknesses 
of b = 1 nm and b = 2 nm in Figs. 3(a) and 3(b), respectively. The solid line represents 
the exact calculation, which was obtained using Eq. (J3J). The dotted line is the result of 
our analytical simple model described by Eq. (|14p . Note that the extremal point for the 
NM Fermi surface coincides with the maximum transmission point of the tunnel barrier at 
k|| = 0. When the NM layer is thin, there is some discrepancy between the exact result and 
our analytical model, but the agreement improves as the NM thickness increases. The overall 
trend of the NM thickness dependence is well depicted by our analytical model. The TMR 
oscillates and goes to zero as the NM thickness increases. It is also shown that the TMR 
dependence on d is affected by the thickness of the tunnel barrier. The peak points of TMR 
do not coincide for different tunnel barrier thicknesses. The decay of TMR as a function of 
d is faster for a thinner tunneling barrier. It is essential to consider the kii dependence of 
the transmission coefficient in the calculation. The effect of the tunnel barrier thickness on 
the NM thickness dependence of TMR is well described by our proposed analytical formula 
Eq. AH]). 



B. Spacer with nonparabolic dispersion relation 

Second, to investigate the case of the multiple extremal spanning vectors in the NM layer, 
we assumed the following effective dispersion relation for the NM layer: 
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(15) 



where the constant a is set to a = 0.68. Except for the NM layer and b = 1.5 nm, the 
dispersion relations and parameters for the FM and tunnel barrier are the same as in the 
previous case. The effective mass and k z in the NM layer were determined using Eq. ffl5|) . 
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In Fig. [1(a), the cross section of the Fermi surface for bulk NM is plotted as a function of 
kn. There are two kinds of extremal spanning vectors : one at fcii = and the other at 
fey = 0.68/cp, where kp = ^2moEp/h is the magnitude of the Fermi wave vector for the 
FM with majority spin. The extremal spanning vector at kn = (kn = 0.68kp) is shorter 
(longer) and can give rise to a long (short) period of the TMR oscillation, which is similar 
to the Cu(OOl) case. The TMR as a function of d is shown in Fig. 11(b). The solid line is the 
exact calculation and the dotted line is based on Eq. f JT4|) . The agreement is fairly good, and 
only oscillation with a long period was observed. For the analytical model calculation, we 
included only the extremal spanning vector at kn = and ignored the contribution from the 
extremal point at k\\ = O.Q8kp. This is because the tunneling probability decreases rapidly 
away from k\\ = and the spin asymmetry of the transmission coefficient (|t^| 2 — \t^\ 2 ) is 
very small at fcii = 0.68/cp. The observed behavior clearly shows that even though there 
are two possible periods of the TMR oscillation from the Fermi surface of the NM, only the 
period with significant spin asymmetry of the transmission coefficient would survive. The 
case of multiple extremal spanning vectors in the NM layer is also described well by our 
analytical formula. 



C. Tunnel barrier with nonparabolic dispersion relation 

Third, we studied the case that the point for the maximum tunneling probability does 
not coincide with the position of the extremal spanning vector of the NM Fermi surface in 
k|j space. We assumed that the dispersion relation of the tunnel barrier is given by 

_ \ h\k 2 x - a 2 k 2 F f h\k 2 y - a 2 k 2 F ) 
Ami ^ m o 

where Vb is the bottom energy of the tunnel barrier, V\\ is an energy parameter, and kp is 
the magnitude of the Fermi wave vector for the NM material. With this tunnel barrier, the 
tunneling probability is highest at kn = (akp, akp). For the calculation, we used b = 1.5 nm, 
Vb = 6 eV, V\\ = 0.7 eV, and a = 0.566. The other materials are assumed to be the same as 
in Sec. III-A. The dispersion relation of the FM material is given by e(k) = h 2 k 2 /2m^ + V a 
with = 0, Vj. = A, m-f = m , and = moEp/(EF — A). The dispersion relation of the NM 
material is the same as that of the FM material for majority spin. Note that the extremal 
spanning vector of the NM Fermi surface is located at ky =0. The parameters used in the 
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calculation are ep = 4 eV and A = 2.5 eV. The TMR dependence on the thickness of the 
NM layer is displayed in Fig. 5. The solid line is the exact result and the dotted is the result 
of our analytical formula. As expected, the TMR decays faster and it is almost negligible 
when the NM layer is thicker than about 2 nm. The inset is the transmission coefficient as 
a function of kii along the [110] direction for the FM/I/NM system. It has the maximum 
value around fcno = 0.8kp, at which we have Vky (x) = 0. On the other hand, the spanning 
vector (q) of the NM Fermi surface has the longest at kn = 0. The maximum point of the 
tunneling probability does not coincide with the extremum of the NM Fermi surface, and 
the extremal point is determined from Vk, [— b x + i{qd + <p)] = 0. The extremal point is 
(^-dependent and a complex number. It is not easy to calculate the (^-dependent extremal 
point exactly and we used the following approximation. At the Fermi level, k z in Eq. ( TTSj) 
is an imaginary number [k z = in). We used a parabolic function of k which was expanded 
in a Taylor series around kii = (akp,akp). The calculation became much simpler, and the 
extremal point and the corresponding vector q a in Eq. ( fl3l) were obtained immediately as 
functions of d. In the range of thin NM layers with significant TMR, q a is shorter than 
2fcp, the extremal spanning vector of the NM Fermi surface. Thus, the oscillation period 
in this region is longer than what is expected from the NM Fermi surface. Whenever the 
extremal spanning vectors of the NM Fermi surface are significantly away from the point of 
the maximum transmission coefficient in the kn space, we expect rapid decay of TMR as a 
function of the NM thickness. 

D. Free-electron model 

Finally, we investigated the free-electron model case, where all the effective masses are 
simply the bare electron mass mo- The NM band is assumed to be the same as the majority- 
spin band of FM, which is commonly adopted in the theoretical calculations. The parameters 
used in the calculation are Bp = 4 eV, V[ = 6 eV, A = 3.5 eV, and b = 1 nm. A plot of 
the TMR dependence on NM thickness is displayed by a dotted line in Fig. 6. In this case, 
the TMR reaches a finite value when the NM thickness becomes infinite. At first sight, this 
seems to contradict our analytical formula Eq. ( fl~4l) given in Sec. II. The reason for finite 



TMR at infinite d can be explained as follows. We obtained k z = y 2m (eF — A)/h 2 — fc| 
from the dispersion relation at the Fermi level. When the magnetizations of the FM layers 
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are antiparallel, electrons with the majority spin in the FM(L) layer do not penetrate to the 
FM(R) layer for h\ > ^/2mo(ep — A)/K because k z becomes imaginary in the FM(R) layer. 
The electrons with the majority spin for k\\ > ^/2mo(ep — A)/h in the FM(L) layer are 
totally reflected at the NM/FM(R) interface due to the potential step and do not contribute 
to conduction. Thus, Gap is underestimated and gives rise to a finite AG and consequently 
a finite TMR, even when the NM layer is infinite, as long as the transport is coherent. 
This is more pronounced for a larger A, lower barrier height, and thinner tunneling barrier 
because a larger portion of electrons with the Fermi energy will be completely reflected at the 
NM/FM(R) interface. When deriving Eq. ( |T3|) in Sec. II, we considered multiple reflection 
inside the NM layer to calculate the transmission, and we assumed that once the electrons 
tunnel through the tunnel barrier, most of them flow to the FM(R) layer after multiple 
reflections. However, in the AP magnetizations of FM layers of the free-electron model, 
the electrons with majority spin for k\\ > y^2m (eF — A)/h in the FM(L) layer cannot 
penetrate into the FM(R) layer after tunneling because of the potential step, resulting in a 
finite TMR for the infinite NM layer. To clarify this point, we calculated the TMR for the 
finite FM(R) layer; namely, the FM(L)/I/NM/FM(R)/NM junctions. The FM(R) layer is 
1 nm thick, and the result is shown by the solid line in Fig. 6. The FM(R) layer behaves 
like a potential barrier for electrons with the majority spin and k\\ > ^/2m (eF — A)/h in 
the FM(L) layer. However, the TMR goes to zero as the NM thickness increases. When 
the FM(R) layer is not too thick, the tunneling probability through this potential barrier 
is larger than that through the tunnel barrier (I), and most of electrons with majority spin 
and k\\ > ^/2mo(ei? — A)/h in the FM(L) layer eventually flow to the NM(R) layer once 
they tunnel through the insulating barrier. Even for the infinite FM(R) layer, the TMR for 
the infinite NM layer becomes negligible as A increases, the thickness of the tunnel barrier 
increases, and the barrier height increases. This is because only a small portion of electrons 
with the Fermi energy will be completely reflected at the NM/FM(R) interface. Also, if the 
minority spin band of the FM layer is set to be the same as the NM band, the TMR goes 
to zero as the NM thickness increases because there is no total reflection at the NM/FM(R) 
interface. The finite TMR for the infinite NM layer is possible when a significant portion of 
electrons with the Fermi energy have less transmission probability into the FM(R) layer than 
that through the tunnel barrier. Complete reflection of electrons at the interface between 
two metals may occur for some given kii at the Fermi level due to the mismatch of the 
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electronic states. However, even slight scattering would lead to the penetration of electrons 
into the FM(R) layer. In this sense, we do not expect that a large portion of electrons 
with Fermi energy are reflected completely, in reality, at the NM/FM interface. Thus, it is 
unrealistic to expect a finite TMR as the NM thickness increases in experiments. 

IV. CONCLUSION 

We calculated the TMR of FM/I/NM/FM tunnel junctions. The TMR was calculated as 
a function of NM thickness using the Landaur-Biittiker formula. Multiple band structures 
were included and an analytical form describing the TMR was obtained. Conductance 
was calculated from the summation of the transmission over ky. The transmission was 
obtained by considering multiple reflections between the I/NM and NM/FM interfaces. 
The summation over kn was carried out analytically. The contribution was mainly from 
the extremal point in ky-space that was determined from the combination of the NM Fermi 
surface and the ky -dependence of the transmission coefficient of the FM/I/NM junction. The 
TMR was expressed with the transmission coefficient of the FM/I/NM junction, reflection 
amplitudes at the NM/I and NM/FM interfaces at the Fermi level, and the extremal wave 
vector. Many oscillation periods can be inferred from the shape of the NM Fermi surface, 
but they can be observed only when the corresponding spin asymmetry of the transmission 
coefficient is significant. We suggest that only few oscillation periods are likely to be observed 
in real experiments. When the NM spacer is thin, our proposed model indicates that the 
decay of the TMR was slower than the inverse of the space thickness for coherent transport. 

Numerical calculations were performed to investigate the accuracy of the proposed for- 
mula. An envelope-function theory was adopted, and our model was compared to the exact 
result. We showed that the results of the proposed formula are in good agreement with the 
exact calculations. The TMR dependence on the thickness of the NM spacer was affected by 
the tunnel barrier thickness, which was well described by our formula. The numerical calcu- 
lation was extended to the case with multiple extremal spanning vectors in the Fermi surface 
of the NM spacer. Our proposed formula is in good agreement with the exact result, and 
only the oscillation period with significant spin asymmetry of the transmission coefficient 
was observed as predicted using our formula. When the tunneling probability associated 
with each extremal spanning vectors of the NM Fermi surface is low, the TMR decays very 
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fast as the NM thickness increases. A free-electron case was also considered. The NM band 
was assumed to be the same as the majority spin band of the FM layer, and spin-splitting 
in the FM layer was assumed to be rather large. The TMR approached a finite value as the 
NM thickness increased. This was because a large portion of electrons with the majority 
spin in the left FM layer were reflected completely at the right NM/FM interface. When 
the semi-infinite FM layer on the right side was replaced by a finite layer, the TMR decayed 
to zero as the NM thickness increased. As long as the transmission into the FM layer was 
higher than that through the tunnel barrier, the TMR became zero as the thickness of the 
NM spacer increased. We suggest that finite TMR for the infinite NM spacer is unrealistic 
in real experiments. 

Acknowledgments 

This work was supported by the Basic Science Research Program through the National 
Research Foundation of Korea (NRF), funded by the Ministry of Education, Science and 
Technology (Grants No. 2011-0005251 and No. 2012-0002984). 



* Electronic address: chan@i nha.ac.krl 

1 J. S. Moodera, L. R. Kinder, T. M. Wong, and R. Meservey, Phys. Rev. Lett. 74, 3273 (1995). 

2 S. Yuasa and D. D. Djayaprawira, J. Phys. D: Appl. Phys. 40 R337 (2007). 

3 A. Vedyayev, N. Ryzhanova, C. Lacroix, L. Giacomoni, and B. Dieny, Europhys. Lett. 39, 219 
(1997). 

4 P. LeClair, H. J. M. Swagten, J. T. Kohlhepp, R. J. M. van de Veerdonk, and W. J. M. de 
Jonge, Phys. Rev. Lett. 84, 2933 (2000). 

5 P. LeClair, J. T. Kohlhepp, H. J. M. Swagten, and W. J. M. de Jonge, Phys. Rev. Lett. 86, 
1066 (2001). 

6 S. Zhang and P. M. Levy, Phys. Rev. Lett. 81, 5660 (1998). 

7 S. Yuasa, T. Nagahama, and Y. Suzuki, Science 297, 234 (2002). 

8 H. Itoh, J. Inoue, A. Umerski, and J. Mathon, Phys. Rev. B 68, 174421 (2003). 

9 H. Itoh, J. Inoue, A. Umerski, and J. Mathon, J. Magn. Magn. Mater. 272-276, el467 (2004). 



15 



A. A. Shokri and A. Saffarzadeh, J. Phys.:Condens. Matter 16, 4455 (2004). 

Z. M. Zeng, X. F. Han, W. S. Zhan, Y. Wang, Z. Zhang, and S. Zhang, Phys. Rev. B 72, 054419 
(2005). 

J. Yang, J. Wang, Z. M. Zheng, D. Y. Xing, and C. R. Chang, Phys. Rev. B 71, 214434 (2005). 
Z. P. Niu, Z. B. Feng, J. Yang, and D. Y. Xing, Phys. Rev. B 73, 014432 (2006). 
X. Feng, O. Bengone, M. Alouani, I. Rungger, and S. Sanvito, Phys. Rev. B 79, 214432 (2009). 
S.-P. Chen and C. R. Chang, IEEE Trans. Magn. 45, 2410 (2009). 

G. Autes, J. Mathon, and A. Umerski, Phys. Rev. B 80, 024415 (2009). 
S.-P. Chen, J. Appl. Phys. 107, 09C716 (2010). 

S.-P. Chen, Thin Solid Films 519, 8215 (2011). 

B. C. Lee, J. Appl. Phys. 107, 09C708 (2010). 

H. Bruus and K. Flensberg, Many-Body Quantum Theory in Condensed Matter Physics (Oxford 
University Press, New York, 2004). 

B. Lee and Y.-C. Chang, Phys. Rev. B 52, 3499 (1995). 

B. C. Lee and Y.-C. Chang, Phys. Rev. B 62, 3888 (2000). 

M. D. Stiles, J. Magn. Magn. Mater. 200, 322 (1999) and references therein. 
L. Tsetseris, B. C. Lee, and Y.-C. Chang, Phys. Rev. B 55 11586 (1997). 

C. -Y. You, C. H. Sowers, A. Inomata, J. S. Jiang, S. D. Bader, and D. D. Koelling, J. Appl. 
Phys. 85, 5889 (1999). 

T. Nagahama, S. Yuasa, E. Tamura, and Y. Suzuki, Phys. Rev. Lett. 95, 086602 (2005). 

R. Matsumoto, A. Fukushima, K. Yakushiji, S. Nishioka, T. Nagahama, T. Katayama, Y. 

Suzuki, K. Ando, and S. Yuasa, Phys. Rev. B 79, 174436 (2009). 



16 



(a) 



A + — 






— B 


A 






<— B 


FM(L) 


I 


NM 


FM(R) 



< >< > 

b d 









(c) 




^B* + 


A L ~* — 


I 








^B* 


FM(L) 


NM 




NM 


FM(R) 



< > 

b 



FIG. 1: (a) Schematic diagram of a magnetic tunnel junction (MTJ) with a nonmagnetic (NM) 
layer inserted between the insulating (I) tunnel barrier and the right ferromagnetic layer [FM(R)]. 
d is the thickness of the NM layer. A + , A - , B + , and are coefficient vectors [see Eq. (pQ)]. The 
transmission in the MTJ can be expressed with the reflection and transmission amplitudes of the 
separated interfaces shown in (b) and (c). 
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FIG. 3: TMR as a function of d for the effective-mass band. The barrier thicknesses are (a) b = 1 
and (b) b = 2 nm. The parameters used in the calculation are ep = 4 eV, Vj = 6 eV, and A = 2.5 
eV. The effective mass of electrons with minority spin in the FM material is m| = 77?,o£,f/(£f — A) 
and other effective masses are the bare electron mass mo. The solid line is the exact result, and 
the dotted line is based on the proposed analytical formula. 
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FIG. 4: (a) Cross section of the Fermi surface for the NM spacer with the dispersion relation given 
in Eq. (|15p . kp = \J2niQEF /h is the Fermi wave vector for the majority spin of the FM layer. It is 
similar to the Cu(001) case, (b) TMR as a function of the NM layer thickness d. Except for the 
NM spacer, the other parameters are the same as those used in the previous case. The solid line 
is the exact result, and the dotted line is based on the proposed analytical formula. 
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FIG. 5: TMR as a function of the NM layer thickness d when the point for the maximum trans- 
mission coefficient of the tunnel barrier does not coincide with the extremal point of the NM Fermi 
surface in ky space. The dispersion relation of the tunnel barrier is given in Eq. (|16p . and the FM 
and NM materials are assumed to be the same as in Fig. 3. The inset is the transmission coefficient 
as a function of kn along the [110] orientation for the FM/I/NM system. 
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FIG. 6: TMR as a function of the NM layer thickness d for the free-electron band. The parameters 
used in the calculation are ef = 4 eV, Vj = 6 eV, A = 3.5 eV, and 6 = 1 nm. The solid line is for 
the 1 nm FM layer and the dotted line is for the infinite FM layer. 
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